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ABSTRACT 

The volume coefficient K (=incompressibility of the nuclear matter), the Coulomb 
coefficient K^, and the volume-symmetry coefficient Kyg of the nucleus incom¬ 
pressibility are studied in the framework of the relativistic mean field theory, 
with aid of the scaling model. It is found that K = 300 ± 50MeV is neces¬ 
sary to account for the empirical values of Ky, Kc, and Ky^, simultaneously. 
The result is independent on the detail descriptions of the potential of the a- 
meson self-interaction and is almost independent of the strength of the w-meson 
self-interaction. 


One way to determine the incompressibility K of nuclear matter from the 
giant monopole resonance (GMR) data is using the leptodermous expansion[l] 
of nucleus incompressibility K{A, Z) as follows. 

K{A,Z) = K + KsfA-^/^ + K^sI^ + K^Z^A-^/^ + --- ; I=l-2Z/A,(\) 

where the coefficients, K, Kgf, and are volume coefficient (incompress¬ 

ibility of nuclear matter), surface term coefficient, volume-symmetry coefficient 
and Coulomb coefficient, respectively. We have omitted higher terms in eq. (1). 
Although there is uncertainty in the determination of these coefficients by us¬ 
ing the present data, Pearson [2] pointed out that there is a strong correlation 
between K and Kg- (See table I.) Similar observations are done by Shlomo and 
Youngblood [3]. 


Table I 


According to this context, Rudaz et al. [4] studied the relation between 
incompressibility and the skewness coefficient by using the generalized version 
of the relativistic Hartree approximation [5]. The compressional and the surface 
properties are studied by Von-Eiff et al. [6-8] in the framework of the relativistic 
mean field approximation of the a-uj-p model with the nonlinear cr terms. They 
found that low incompressibility {K « 200MeV) and a large effective nucleon 
mass M* at the normal density (0.70 < M*/M < 0.75) are favorable for the 
nuclear surface properties [8]. On the other hand, using the same model, Bodmer 
and Price [9] found that the experimental spin-orbit splitting in light nuclei 
supports M* « 0.60M. The result of the generator coordinate calculations for 
breathing-mode GMR by Stoitsov, Ring and Sharma [10] suggests K w 300MeV. 

In previous papers [11,12], we have studied the relation between K and Kc 



in detail, using the relativistic mean field theory with the nonlinear a terms 
[13] and the one with the nonlinear a and w terms [14]. We found that, under 
the assumption of the scaling model [1], K = 300 ± 50MeV is favorable to 
account for K, Kc and Ky^, simultaneously. It seems that this conclusion is 
not drastically changed in the use of the relativistic mean field theory and the 
scaling model. In this paper, we examine the conclusion in more general way, 
in which the result does not depend on the detail descriptions of the cr-meson 
self-interaction. The reason why we restrict our discussions to K, Kc and Kyg is 
that the general discussions, which are independent of the detail of the model ( 
e.g., types of the interactions, values of the parameters in the Lagrangian, etc.) 
are possible to a considerable extent, since, as is shown below, these quantities 
are almost analytically estimated by using the result for the nuclear matter, if 
we assume the scaling model [1]. 

We use the relativistic mean field theory based on the a-LO-p model with the 
nonlinear cr terms. The Lagrangian density consists of four fields, the nucleon ip, 
the scalar cr-meson p, the vector w-meson V^, and the vector-isovector p meson 
b^, i.e., 

Lnctlop = - M)'tp 

+ 

+gs'ip'tp(t> - - gP'4’1 ; 

gp^^ X t)j/, (2) 

where Wg, rriy, nip, ps, Pv and Pp are ct- meson mass, w-meson mass, p-meson 
mass, cr-nucleon coupling, w-nucleon coupling, and p-nucleon coupling, respec¬ 
tively. The U{(p) is a nonlinear self-interaction potential of a meson field </>. For 
example, in ref. [11], we have used the quartic-cubic terms of (p as in ref. [13], 



( 3 ) 


i.e., 

[/(</>) = ^ 

where b and c are the constant parameters which are determined phenomeno¬ 
logically. However, in this paper, we do not give an explicit expression of U{(j)) 
and discuss the problem in more general way, without any assumption of U{(f>). 
In the scaling model [1], K, Kc and Kyg in eq. (1), are given by 


K = 9pt 
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where p, po, Ei, and q^i are the baryon density, the normal baryon density, the 
binding energy per nucleon, and the electric charge of proton, respectively, and 
Ro = [3/(47rpo)]^/^ 

T^' _ o I 


dJ , 


L — 3po lp=P0! Ksym — 9Pq ^^2 [P=Po • 


dp^ 

'dp-^ 


r 1 
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(7) 

( 8 ) 


The quantity such as K' is sometimes called ’’skewness”. 

In the mean field theory with the Lagrangian (2), L and Ksym are given by 

[7] 


L = 


3p ^2 


8M2 


p ( 2kFkp kpEp 

2\^ ^ 


(P = Po), 
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and 


3 ^f2k'^ , 2fcffc" Akpk'pE*; , 2epE*;^ k^E*;' 

sym 1 7^^ “r ^ Tp ^2 ' Z7'*2 

z \ h/p h/p h/p H/p h/p 

where kF{= [37r2p/2]^/^) is the Fermi momentum, and Cp = 

, dkF kF d'^kF —2kF 

^ dp ip' ^ dp'^ 9p2 


j {p = Po), (10) 

PpM/mp, 


( 11 ) 



E*p = Jkl + M*^, = 


d^E*p 


M* in eq. (12) is the effective nucleon mass. Furthermore, at p = po, Ep and 
Ep' are related to K and K' in the following relations, respectively [7,11]. 


p- _ JL _ ^ 

^ 9po M2 ’ 


E*j^' = 


K + K' 


At p = Po, Cv and Cp are also related to M* as follows. 


= —\^M -ai- ^kl + M*2 j , 


[13] and 


QE*p)' 


[15] where oi and 04 are the binding energy and the symmetry energy at p = po, 
respectively. We remark that eqs. (5), ( 6 ), and (9)~(16) have no explicit 
dependence on U{(j)). 

From eq. (5), K' are determined, if po) K and Kc are given. Therefore, from 
the eqs. (5), ( 6 ), and (9)~(16), it is seen that is determined, if po, M, oi, 04 , 
K, Kc, and M* are given, without giving the detail descriptions for [/(</>). Using 
these equations, we calculate Kyg. In the calculations, we put po = 0.16fm“^, 
M = 939MeV, oi = 15.75MeV and 04 = 30.0MeV. For K and Kc, we use the 
values in table I. We assume that M* = 0.5M ^ 0.94M, the phenomenologically 
acceptable values. (We remark the upper bound for M* is gotten, if we put 
Cl, = 0 in eq. (15). ) In fig. 1 (solid lines), we show Kyg as a function of M* for 
two sets of K and Kc in table I. In the figures, Kyg decreases as M* increases. 


This is understood as follows. As M* increases, L and Ksym decreases, because 
of large Ep in the denominators in eqs. (9) and (10). Small Ksym makes Kyg 
smaller (more negative) and, on the contrary, small L makes Kyg larger (less 



negative), since {9K'/K + 6) > 0 in eq. (6). In the cases of fig. 1, the effect of 
the small Ksym overcomes that of the small L, therefore, Kyg decreases. In the 

case of K = 300MeV and = -3.990MeV, = -289 -98MeV. These 

values are in good agreement with the corresponding empirical values in table 
I. We also remark that the uncertainty of Kyg in changing M* is comparable 
to the empirical error bar of Kyg in table I. In the case of K = 250MeV and 
Kc = —0.7065MeV, = 50 410MeV. These values are somewhat larger 

than the empirical ones. In table 11(a), we summarize the range of the calculated 
Kys for each set of K and Kc in table I. Comparing the table I and table 11(a), 
we see that K = 300 ± 50MeV is necessary to account for K, Kc and Kyg, 
simultaneously. 


Fig. 1(a),(b). Table 11(a),(b) 


We remark the following three points. 

(1) The results are independent of the form of U{(j)), since eqs. (5), (6), and 
(9)~(16) are required for any type of U{<p). 

(2) The question, whether there are coupling parameters, which reproduce 
the set of K and Kc in table I, or not, is still open, and the answer for the 
question depends on the detail descriptions of U{(1)). For example, if we use the 
quartic-cubic potential (3), we could not find the coupling parameters (i.e., gs, 
gv, b, and c ), which reproduce K = 200MeV and Kc = 2.577MeV [11]. (We 
remark that K and Kc are independent on gp in the mean field theory. ) Also, 
using eq. (3), we get the parameter set for K = 300MeV and Kc = —3.990MeV, 
only in the case of M* = 0.83M [11]. In those cases, parameter set is uniquely 
determined or could not be found, since the number of the parameters is not 



larger than the number of the inputs, i.e., K, Kc, and two conditions for the 
saturation. If the higher terms of (p are added to (3), the wider range of M* 
may be available. However, K = 300 ± SOMeV is necessary to reproduce the 
empirical values of K, Kc, and Kyg simultaneously, for any type of U{(j)). 

(3) The results do not have a strong-dependence on pg, ai and 04 , since 
the calculated Kyg is much more sensitive to the ratio K'jK than to those 
quantities. 

Next we add the following attractive term of vector self-interaction (VSI) 
[14,12] to the Lagrangian (2). 

Lvsi = (17) 

, where T is a positive constant which determine the strength of VSI. Although, 
according to this modification, the formalism for calculating Kyg, which is de¬ 
scribed above, is slightly modified, it is still possible to study M*-Kys relation 
without giving the detail descriptions of U{<p), as in the case of no VSI, i.e., 
V = 0. In fig. 1, we also show Kyg as a function of M*, in the cases of j/ = 1 
and y = 5, where y = gppoY/my [12]. In those figures, it is seen that the 
VSI makes Kyg smaller in the small M* region. Kyg is almost independent 
on y in the large M* region, where the value of M* dominates the features 
of the equations of state. As a result, the uncertainty of Kyg in changing 
M* becomes smaller. In the case of AT = 300MeV and Kc = —3.990MeV, 

Ky, = -298 -282(-303 -289)MeV with y = 1.0(5.0). In the case of 

K = 250MeV and Kc = -0.7065MeV, Ky^ = 50 ~ 147(49 ~ 86 )MeV with 
y = 1.0(5.0). In each case, Kyg with y = 5.0 is not much different from that 
with y = 1.0. Kys is hardly changed, if y increases much more. In table 11(b), 
we summarize the range of the calculated Kyg in the cases of j/ = 5.0. Com¬ 
paring this table with table I, we see that K = 300 ± 50MeV is necessary to 



account for the empirical values of K, Kc, and K^s at the same time, as is in 
the case of no VSI. Although the VSI makes K^s smaller (more negative), the 
change of Kyg is comparable to the magnitude of the empirical error bars at 
the empirical available value of K{^ 300MeV), and, therefore, the conclusion is 
hardly changed. The result is also independent on the detail of U{(/)) as in the 
case of no VSI. 

In summary, we have studied K, K^,, and Kyg by using the relativistic mean 
field theories with the nonlinear a term and with the nonlinear a and to terms, 
with aid of the scaling model. It is found that, in both cases, K = 300±50MeV 
is necessary to account for the empirical values of K, Kc, and at the 
same time. The result is independent on the detail descriptions of U{(j)) and is 
almost independent on the strength of VSI. It seems that this conclusion is not 
drastically changed, if we use any type of the relativistic mean-field theory and 
the scaling model, since the calculated Kyg is most sensitive to the ratio K'jK, 
which is adjusted to the empirical values. 
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Table and Figure Captions 


Table I 

The sets of the empirical values of K, Kc and Kyg in the table 3 in ref. [2], 
(According to the conclusion in ref. [2], we only show the data in the cases of 
K = 150 ~ 350MeV.) All quantities in the table are shown in MeV. 

Table II 

Range of the calculated using the sets of K and K^. in table I as inputs, (a) 
The result in the the mean-field theory only with a meson self-interaction, (b) 
The result in the mean-field theory with a and lo mesons self-interactions with 
y = 5.0. : In each table, ’’upper bound”, ’’mean value”, and ’’lower bound” 
mean that the results are obtained by using the upper bound, the mean value, 
and the lower bound of Kc in table I, respectively. All quantities in the table 
are shown in MeV. 

Fig. 1 Kys as a function of M*. (a) The cases of K = 300MeV and 

Kc = -3.990MeV. (b) The cases of AT = 250MeV and Kc = -0.7065MeV: In 
each figure, the solid line is the result in the the mean-field theory only with a 
meson self-interaction, and the dotted and the dashed lines are the results in 
the mean-field theory with cr and lo mesons self-interactions with j/ = 1.0 and 
y = 5.0, respectively. 




Set 1 

Set 2 

Set 3 

Set 4 

Set 5 

K 

150.0 

200.0 

250.0 

300.0 

350.0 

Kc 

5.861 ±2.06 

2.577 ±2.06 

-0.7065 ±2.06 

-3.990 ±2.06 

-7.274 ±2.06 

K 

66.83 ± 101 

-46.94 ± 101 

-160.7 ± 101 

-274.5 ± 101 

-388.3 ± 101 


Table I 




upper bound 

mean value 

lower bound 

K 

= 150 

940~1738 

728~1426 

517~1113 

K 

= 200 

601~1230 

389~ 918 

178- 605 

K 

= 250 

262~ 722 

50~ 410 

-161- 97 

K 

= 300 

-77~ 214 

-289~ -98 

-506—411 

K 

= 350 

-418~-294 

-651—606 

-921—840 

Table 11(a) 




upper bound 

mean value 

lower bound 

K 

= 150 

939-1103 

728- 863 

516- 623 

K 

= 200 

600- 714 

389- 474 

177- 235 

K 

= 250 

261- 325 

49— 86 

-165—154 

K 

= 300 

-80- -63 

-303—289 

-542—500 

K 

= 350 

-451—416 

-691—628 

-931—839 


Table 11(b) 











